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We consider the ground-state entanglement in highly connected many-body systems, consisting of 
harmonic oscillators and spin-1/2 systems. Varying their degree of connectivity, we investigate the 
interplay between the enhancement of entanglement, due to connections, and its frustration, due to 
monogamy constraints. Remarkably, we see that in many situations the degree of entanglement in a 
highly connected system is essentially of the same order as in a low connected one. We also identify 
instances in which the entanglement decreases as the degree of connectivity increases. 
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I. INTRODUCTION 

Entanglement theory has experienced an impressive 
development in the last decade, mainly due to the key 
role quantum correlations play in quantum information 
science. The novel concepts and mathematical methods 
developed in this new research area are beginning to re- 
veal their usefulness also in different contexts. A striking 
example of this tendency is given by the physics of quan- 
tum many-body systems. As an instance, the analysis of 
the role played by entanglement in quantum phase tran- 
sitions allowed for a deeper understanding of this purely 
quantum phenomenon [3, 0, In this scenario, entan- 
glement theory is also giving a fundamental contribution 
in the development of new methods capable of simulating 
efficiently strongly interacting systems @, H, @] . 

Clearly, the correlations between different parts of a 
many-body system originates by their mutual interac- 
tion. In this sense it is natural to expect that the ground 
state of a strongly interacting and connected quantum 
system will exhibit a high degree of entanglement. How- 
ever, this intuition has to be taken cautiously, since the 
shareability properties of quantum correlations are espe- 
cially non trivial and without classical analogue. One 
of the main differences between classical and quantum 
correlations is the so called monogamy of the latter Q . 
In the classical scenario, the fact that two systems share 
some correlations does not prevent them from being cor- 
related with a third party. On the contrary, two maxi- 
mally entangled quantum systems can share no correla- 
tion at all with a third one. More generally, quantum 
correlations are not infinitely sharable, and the more the 
entanglement the less the number of systems with which 
it can be shared. 

Consider two similar Hamiltonians consisting of the 
same interacting terms between pair of particles, the 
only difference being the degree of connectivity. One 
of them, for instance, has only nearest-neighbor interac- 
tions, while the second has also ncxt-to- nearest-neighbor 
interactions. Let us focus on the ground-state entangle- 
ment between two halves of the system. Naively, the 
more connected hamiltonian is expected to have a larger 
entanglement, since there are more bonds connecting the 
two halves. However, in the more connected system, each 
particle has to share the quantum correlations with a 



larger number of particles, so the connecting bonds may 
give a smaller amount of entanglement. Therefore, it is 
unclear which geometry leads to a larger ground-state 
entanglement. 

In this work we analyze the interplay between the 
enhancement of the ground-state entanglement due to 
connections and its suppression due to monogamy con- 
straints. We consider spin-1/2 and infinite dimensional 
(harmonic oscillators) systems of various geometries with 
two-body interactions and focus on the bipartite entan- 
glement between two halves of the system. Remarkably, 
we sec that in many situations the degree of entanglement 
in a highly connected system is essentially of the same or- 
der as in the case of a low connected one. Actually, we 
can even individuate systems for which the entanglement 
decreases as the degree of connectivity increases. 

Before proceeding, let us mention that there exist some 
works studying how the monogamy of entanglement af- 
fects the ground state properties of Hamiltonians with 
nearest- neighbor interactions, see for example [8|, [9(. The 
ground-state entanglement of a highly symmetric and 
connected system, the so-called Lipkin-Meshkov-Glick 
model, was also computed in [Toj . 

The paper is organized as follows. In the next Sec- 
tion we will introduce the systems we are going to an- 
alyze, in particular recalling known results about en- 
tanglement calculations in the different cases considered. 
Then, in Section [ill! wc will show in details how the en- 
tanglement depends on the connectivity of the systems 
themselves. In order to individuate a general behavior 
we will consider a variety of different connectivity con- 
ditions. Specifically, we will analyze both regular and 
random graph configurations, paying particular atten- 
tion to the case of bipartite graphs. The relationship 
between our findings and related results in the context of 
cntanglcment-area laws will be outlined in Sec lIVI Even 
if a way to quantify the action of monogamy in a multi- 
partite setting is still lacking, we will see in Sec|V] that 
the analysis of bipartite monogamy inequalities can shed 
some light onto our findings. We will close the paper 
in Sec I VII by discussing possible implications of our find- 
ings, in particular in the context of classical simulations 
of quantum systems. 
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II. SPIN AND BOSONIC MODELS 

As said, we consider two paradigmatic systems, namely 
interacting spin-1/2 and bosonic particles. Concerning 
the former, we study a system of n spin-1/2 particles 
under the XX Hamiltonian, 

H = Y^u i [oi®oi+<Ti®o 3 v ], (i) 

where a\ (k = x, y, z) denote the Pauli matrices referred 
to the z-th particle. The coupling Uj will be set different 
from zero when the i — j couple directly interacts, i.e. in 
dependence on the topology and connectivity of the sys- 
tem. The actual value of the nonzero t%j will be chosen 
randomly, in order to have averaged properties and avoid 
the dependence of our results on the details of the interac- 
tion. Interactions of the type ([1]) may model highly con- 
nected physical systems, such as quantum spin glasses 
PH . The entanglement between two parts of the sys- 
tem will be measured by the entropy of entanglement E, 
namely the von Neumann entropy S(p) = —Tr[p log 2 p] 
of one of the reduced subsystems [l2j • 

Concerning the bosonic case, we consider systems con- 
sisting of n harmonic oscillators with quadratic cou- 
pling. Such systems may model discrete versions of 
Klein-Gordon fields, or vibrational modes in crystal lat- 
tices, ion traps and nanomechanical oscillators. We de- 
fine the vector R of quadrature operators by Rj = Xj 

and Rn+j = Pj (1 < j 5= n), where Xj and Pj are the 
position and linear momentum operator respectively. For 
simplicity, we consider only a coupling via the different 
position operators, in which case the Hamiltonian is of 
the form 

where 1„ denotes the n x n identity matrix. The po- 
tential matrix V is defined via the harmonic coupling 
between oscillator i and j, namely a(Xi — Xj) 2 /2. For 
each geometry considered in the following, we denote by 
C the nxn adjacency matrix of the corresponding graph, 
with elements Cjj = Cji = 1 if the i-th and j-th oscilla- 
tor are coupled and cy = Cji = otherwise. Then, the 
potential matrix V is given by Vij = —acij (i ^ j) and 
Vu — 1 + a J2j c ij ■ The ground state of the system is 
a Gaussian state characterized by the covariance matrix 
7 = (7x©7p)/2, with lx = V- 1 ' 2 and lp = V 1 ' 2 0. 
We use as entanglement measure the logarithmic nega- 
tivity Ni, between two gener ic g roup A and B. It 
can be shown that Ni is given by [13| 

n 

iV ; = -^log 2 min[l,A J (7 x P 7p P)], (3) 
j=i 

where Aj(M) is the j-th eigenvalue of matrix M. We 
denote by P the nxn diagonal matrix with j-th diag- 
onal entry given by 1 or —1, depending on whether the 



oscillator on position 1 < j < n belongs to group A or 
B, respectively. 



III. ENTANGLEMENT IN DIFFERENT 
TOPOLOGIES 

A. Neighbor coupling 

The first configuration that we consider is a one- 
dimensional (ID) chain of n particles, in which each of 
them can interact with n c of its neighbors on each side. 
Thus, n c is the parameter that characterizes the degree 
of connectivity in this setting (see Fig. [1]). We consider 
a distance-independent interaction, in order to avoid any 
dependence on the particular scaling of the interaction 
strength. In particular, given the ground state, we cal- 
culate the entanglement between the two halves of the 
system (groups A and B) as a function of the number of 
interacting neighbors n c . The typical behavior in the case 




FIG. 1: Closed chain of n = 10 interacting particles in 
the nearest- neighbor (n c = 1, left), next-to-nearest neighbor 
(n c = 2, center) and fully connected configuration (n c = 5, 
right). The entanglement is computed between two halves of 
the system. 

of a XX system, is reported in Fig.[5]for the case of n = 22 
spins. The exact calculation of the ground state was per- 
formed using the SPINPACK package [3. The solid 
line represents the averaged ground-state entanglement, 
where the Hamiltonian parameters Uj between pair of 
particles are randomly chosen in the interval [0, 1], while 
the dashed line gives the largest entanglement obtained. 
We clearly see that the entanglement grows only slightly 
and, in particular, the fully connected chain has a de- 
gree of entanglement comparable to the nearest-neighbor 
coupled chain. Note that, by contrast, the number of 
bonds connecting the two halves of the chain increases 
as n 2 . The same behavior has been observed for differ- 
ent Hamiltonian operators, consisting of other interaction 
terms, and smaller sizes. 

We consider now the same configuration for the case of 
a chain of harmonic oscillators. As said above, the inter- 
actions between the particles simply correspond to oscil- 
lators of coupling constant a. The entanglement between 
the two halves of an open chain consisting of n = 100 os- 
cillators is shown in the inset of Fig. ® where the logarith- 
mic negativity is plotted versus the number of coupled 
neighbors, n c . One clearly see that the entanglement in- 
creases (almost linearly) as far as n c < n/2, whereas for 
higher connected systems the entanglement is frustrated. 
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FIG. 2: For a closed chain of n — 22 spin-1/2 particles with 
XX interaction ([]]), the entropy of entanglement is plotted 
versus the number of connected neighbors n c , averaged over 
100 realizations. The dashed line gives the largest entangle- 
ment obtained. Inset: for an open chain of n = 100 oscillators 
the logarithmic negativity Ni is plotted versus n c . From top 
to bottom the coupling constant a is given by a — 10, 1,0.1. 



The frustration mechanism is indeed stronger than in the 
spin case, the entanglement decreasing at some point as 
the number of connections increases. Notice the quite 
universal behavior of these curves: the position of the 
maximum does not depend on the coupling constant a 
and, as one can expect, the entanglement increases with 
a, for fixed n c . 

Both the examples reported here confirm that the 
monogamy of entanglement plays a predominant role for 
highly connected systems. As said, as the connectivity 
increases, each particle of, e.g., set A becomes as well en- 
tangled with many other particles of the same set. This 
in turn limits, for monogamy reasons, the entanglement 
with the particles of set B. 



our results, since the degeneracy is equal to zero already 
for small values of the connectivity parameter. Actually, 
this degeneracy disappears for values of c p close to n/n t 
(where nt is the number of connections of the complete 
graph), i.e. when we are able to construct a connected 
graph with high probability. In general we see again that 
the entanglement saturates, and the saturation value is 
reached for low value of the connectivity, showing again 
the strong effect of monogamy. 
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FIG. 3: Entropy of entanglement for a randomly coupled sys- 
tem with n = 10 spins (solid), its degeneracy (dashed-dot) 
and its energy (dashed), as a function on c p . We averaged 
over 100 executions. 

In the case of bosonic oscillators, one also retrieves 
the general behavior already highlighted in the previous 
configuration. Namely, the entanglement increases up 
to a certain value of the connectivity parameter [c p in 
the present configuration), after which the monogamy 
constraint starts causing a decrease of the entanglement. 



B. Random coupling 

We have seen above that the details of the entangle- 
ment frustration mechanism depend on the system un- 
der consideration, nevertheless the general behavior is 
the same for a variety of situations. As another example, 
particularly different from the ID chain above, consider 
now a random configuration in which each couple of par- 
ticles i and j is connected with probability c p , which will 
be called connectivity parameter. Then, the connections 
are given by the edge of a graph with random (symmet- 
ric) adjacency matrix. Concerning the case of the XX 
system, the typical behavior is reported in Fig. [3l where 
the entropy of entanglement is plotted as a function of 
the connectivity parameter. We also plot the degener- 
acy of the ground state. Indeed, if the degeneracy was 
non-zero for a significant range of values of c p , the re- 
sults should be carefully considered, since there may be 
other ground states having quite different entanglement 
properties. However, this does not affect the validity of 



C. Bipartite graphs 

We exploit now the configuration in which perhaps the 
effects of monogamy show up more impressively, namely 
the case in which the system can be represented by a 
bipartite graph. The latter is constituted by two sets, 
A and B, for which particles belonging to the same set 
do not directly interact (see Fig. [4]) . Then, a priori, one 
could expect that in such a configuration the effects of 
monogamy should be weaker than in the configurations 
discussed above. We report here the results for both reg- 
ular and random bipartite graphs in the case of harmonic 
oscillators. 

Let us first consider the case of a random bipartite 
graph. As above, we name connectivity parameter c p the 
probability that a generic particle in A interacts with a 
particle in B. For example in the fully connected case 
(c p = 1) each particle in A is coupled to all the particles 
in B. For a fixed coupling constant we look for the op- 
timal c° pt such that the entanglement is maximized. As 
shown in Fig. we have that c° pt ^ 1 in general, de- 
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FIG. 4: Regular bipartite graph of n = 10 interacting parti- 
cles for increasing connectivity. 



pending non-trivially on a. Namely, for large values of a 




The eigenvalues of Vb g arc given by 



1 + 2an c - 2a Vcos[(2j - l)k — 1, (5) 



hence, following Ref. [l3[, the log- negativity between the 
sets A and B can be straightforwardly obtained: 



n/2-l 
fc=0 



log 2 



A* 



A*-* 



(6) 



When n is large, the sum above turns into a Riemann 
series, which gives an explicit expression for the log- 
negativity: 



N, 



(7) 



where we have defined the function 



/(a, n c ) = dx 



log 2 |l + 2an c - 2a^cos[(2j - 



lo: 



g 2 {l + 2an c - 2a^ cos[(2j - 1)(tt - a?)]} 



(8) 



FIG. 5: For a random bipartite graph of harmonic oscillators 
the logarithmic negativity Ni is plotted (normalized for each 
a) versus the connectivity parameter c p , for different coupling 
constants a. The total number of particles is n = 100. 

the maximum entanglement is provided by an Hamilto- 
nian with few connections for each oscillator. Vice-versa, 
for low values of a the completely connected case tends 
to maximize the entanglement. Note again that parti- 
cles belonging to the same set do not directly interact. 
However, these particles become entangled through the 
common interaction with the particles of the other set. 
This, at the same time, limits the amount of entangle- 
ment between the two sets because of monogamy. 

In order to obtain a deeper insight into the behavior 
described above, let us now consider a regular graph. In 
this case, remarkably, it is straightforward to obtain an 
explicit expression of the log-negativity for an arbitrary 
large number of particles. Let us focus on a periodic 
system consisting of n oscillators (labeled from 1 to n) 
and consider two sets composed such that the even os- 
cillators belong to A and the odd ones to B. If each 
oscillator in set A interacts only with 2n c oscillators in 
B (specifically, the 2j-th oscillator interacts with oscil- 
lators labeled 2k + 1 with k = j — n c , . . . ,j + n c ), then 
the interaction matrix V in Eq. ([2]) is given by a circu- 
lant matrix Vbg, whose first row is given by the following 
vector 



(l+2cm c , —a, 0, 



-a, 0,0, 



n c times 



,0,0, -a,..., 0, -a) 

n c timcs 

(4) 



in order to single out the dependence on a and the con- 
nectivity parameter n c . Noticeably, in Eq. ([7]) the depen- 
dence on n factories, allowing to analyze how the con- 
nectivity affects the entanglement by simply looking at 
/(a, n c ) for different coupling regimes (sec Fig.[H]). It can 
be seen that the results arc very similar to the ones of a 
random bipartite graph, indicating that the asymmetries 
of the latter do not play a significant role. As in Fig.O for 
weak coupling the entanglement increases with the con- 
nectivity of the system, meaning that monogamy does 
not play a significant role. On the other hand, for strong 
coupling the entanglement suddenly decreases with n c . 
This is a clear sign of the action of monogamy constraints, 
since for high values of a a high degree of entanglement 
is created already for small n c . 



IV. CONNECTIVITY AND AREA LAW 

In this section we want to elucidate the relation be- 
tween the findings exposed above and the general issue 
concerning entanglement-area laws. A number of the- 
oretical investigations has shown that, for a variety of 
physical systems, the entanglement between two comple- 
mentary regions scales as the area between them |16j |. 
For instance, it is shown in Refs. 0, [TtJ that, for non- 
critical systems with nearest-neighbor coupling, the en- 
tanglement between a distinguished part of a system and 
the rest increases as the number of connections. Regard- 
ing more general interactions, in particular long-range 
ones, the boundary area (given by the connections be- 
tween the two regions of the system) only gives an upper 
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FIG. 6: For a regular bipartite graph of harmonic oscillators 
the function f(a, n c ) in Eq. © is plotted versus the connec- 
tivity parameter n c (arbitrary units, normalized for each a), 
for different coupling constants a (see text for details). 

bound to the entanglement In all these works the 

Hamiltonian of the global system is kept fixed, whereas 
the size of the distinguished region varies. It is then clear 
that the monogamy constraints do not act significantly, 
since the connectivity of the systems remains unchanged 
with the size. On the contrary, in most of our previ- 
ous analysis we kept fixed the size of the distinguished 
region, whereas the connections between particles were 
modified by changing the degree of connectivity of the 
Hamiltonians (see Figs. [21 [3] and [5}. 
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configuration exploited in Scc. lIII A1 In particular, we fo- 
cused on the results corresponding to i) nearest-neighbor 
coupling and ii) the optimal configuration in which the 
number of connections n c is chosen in order to give the 
maximal amount of entanglement. For a closed harmonic 
chain, wc observe that the entanglement remains con- 
stant in the nearest-neighbor case (as we expect from the 
results of Ref . [l3| ) , whereas it increases only logarithmi- 
cally in the optimally connected case. Actually, an exact 
logarithmic increase can be derived in a fully connected 
topology (i.e., n c = n), for which the computation of the 
log-negativity is straightforward following again Ref. . 
There, a lower bound to the log-negativity is given as a 
function of the coupling parameters. Fortunately, such 
bound turns out to be tight in the case of both nearest- 
neighbor and fully connected systems [l9| . In particular, 
for the latter one the log-negativity is given by 

iV, = ilog a [l + 2na] , (9) 

proving the logarithmic increase with n. Remarkably, the 
behavior is similar for a closed spin chain, as can be seen 
in Fig. [Jj Recall that in this case the optimal number of 
connections n° pt is always given by n° pt = n, as pointed 
out in Sec. IIII Al Although our computations are not 
very conclusive, they suggest a sub-linear increase in this 
case too. 

As said, the scenario in Fig.[jjis suitable for a compar- 
ison with the works concerning area laws, in particular 
the ones dealing with half spaces (see, e.g., Ref. [20l]). 
We see that, as expected, for nearest neighbor inter- 
actions an entanglement-area law is satisfied (since we 
are not dealing with critical systems). On the contrary, 
the slight entanglement increase for highly connected sys- 
tems strongly contrasts with the increase of the number 
of bonds linking the two halves of a chain, which scales as 
(n/2) 2 . As recalled, in these cases the boundary area only 
gives an upper bound for the entanglement. Thus, our 
results reveal that the entanglement can actually scale 
sensibly slower than the area in highly connected sys- 
tems. 



V. MONOGAMY INEQUALITIES 



FIG. 7: For an XX spin system the entropy of entanglement E 
between the two halves of the system is plotted as a function of 
the system size n (averaged over 100 realizations). Inset: cor- 
responding graph for the negativity TV; in a closed harmonic 
chain with a = 1 for nearest-neighbor coupling (dashed line) 
and optimal coupling (solid line, see text). 

Following an approach suitable for a comparison with 
the works on area laws, we exploited also the dependence 
of the entanglement on the size of the system. The re- 
sults are reported in Fig. [7J where the main panel refers 
to the case of an XX spin system and the inset to an har- 
monic one. In both cases we considered the non-random 



To quantify how the monogamy of entanglement acts in 
the analyzed scenarios we can refer to monogamy inequal- 
ities. In particular, for the case of harmonic oscillators we 
considered the inequality introduced in Ref. [2l|, where 
the Gaussian tangle r is used as entanglement measure. 
We recall that the latter is defined as the square of the 
negativity, for pure states, and then extended to mixed 
states via a convex roof. In particular, the monogamy 
constraint for n modes can be expressed as 

r , :2 .., i> ^ T Ui (10) 

j=2 
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where r 1:2 - " denotes the tangle between the first mode 
and the remaining ones, whereas t^ 3 is the two-mode 
tangle between the first mode and the j-th one. The dif- 
ference between the lhs and the rhs of Eq. pop is called 
residual entanglement and it indicates the presence of a 
complex structure of the quantum correlations, not sim- 
ply ascribablc to two-mode entanglement. 

The results are reported in Fig. [5] for bosonic oscilla- 
tors in a regular bipartite graph (see Sec. IIII Cj) . We see 
that when the residual entanglement begins to be com- 
parable to the sum of the two-mode tangles, then the 
entanglement between sets A and B starts to be sup- 
pressed. In other words, when a complex entanglement 
structure appears, then particles belonging to the same 
set start to become quantum correlated. As a conse- 
quence of monogamy, they can then share no more a 
high amount of correlations with the particles belonging 
to the other set, and the bipartite entanglement between 
A and B is suppressed. Notice that as the strength of 
the coupling between the oscillators increases the effects 
of monogamy become more and more important, as one 
may expect by general considerations. The additional 
frustration in proximity of the fully connected graph can 
be instead attributed to symmetric reasons. More specif- 
ically, in a fully connected graph the state of the system 
is completely symmetric, which in turn implies a reduc- 
tion of the effective size of the Hilbert space associated 
to the n oscillators. As a consequence, the entanglement 
in the system is frustrated too. 

For spin systems we considered the monogamy inequal- 
ity presented in Refs. 0,113 where r is now the square of 
the concurrence, or tangle, a measure of entanglement 
between qubits. Fig. [5] shows, for the regular bipar- 
tite graph with n c connected neighbors, the quantities 



En 
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and 



.1:2. ..n 



For any value of the connectiv- 



. <7=2 'I 

ity, any spin is maximally entangled with the rest, since 
T i:2...n _ ^ However, as the number of connections 
increases the structure of entanglement becomes highly 
non-trivial. This, as said, limits the bipartite entangle- 
ment between sets A and B. Notice that, due to the 
absence of an on-site interaction in Eq. (fT]), an analogy 
can be seen between the XX system and the harmonic 
one in the strong-coupling regime. 



VI. CONCLUSIONS 

We have analyzed the interplay between ground- 
state entanglement and the connectivity in spin- 1/2 and 
bosonic systems with two-body interactions. We have 
shown that the ground-state entanglement does not nec- 
essarily increase by introducing more interacting terms 
in the Hamiltonian. Actually, for some systems, it does 
decrease with the number of connections. From a more 
applied point of view, the amount of entanglement across 
different bipartitions of a system has been related to the 
complexity of the classical description of the state [23| . 
Although here we just focus on symmetric partitions, 
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FIG. 8: For a regular bipartite graph of n = 90 harmonic 
oscillators the logarithmic negativity Ni is plotted versus the 
number of connections per oscillator n c (solid line, left y- 
axis). The dotted and dot-dashed line are the plot of r 1:2 
and X^j=2 T i 3 respectively (right y-axis). From top to bottom 
the coupling constant is a = 1, 0.1, 0.01. 



it seems reasonable to expect that such partitions have 
the highest entanglement (among the contiguous ones). 
Thus, the previous analysis suggests that the bipartite 
entanglement for many highly connected systems is sim- 
ilar to the one for low-connected Hamiltonian systems, 
where efficient classical algorithms exist, like for example 
Density Matrix Renormalization Group (DMRG). As a 
matter of fact, the latter method has recently been ef- 
ficiently applied to a specific highly connected model, 
in order to analyze quantum phase transitions in spin 
glasses [24j| • The efficiency of DMRG in this scenario 
is not trivial and it is related to the fact that the spe- 
cific system analyzed in Ref. [24[ turns out to be only 
slightly entangled. We have shown here that a large va- 
riety of systems may have such a character, due to the 
fundamental constraint imposed by the monogamy of en- 
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FIG. 9: For a bipartite graph of XX spin systems with n — 14, 
the entropy of entanglement E is plotted versus the number of 
connections n c (solid line, right y-axis). The dashed and dot- 
dashed line are the plot of T 1:2 — n an d ^27=2 T l 3 respectively 
(left y-axis). The average is taken over 100 realizations. 



tanglement. Our results, then, may encourage the search 
for novel classical algorithms able at simulating highly 
connected quantum systems. 
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